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The data of J/ip suppression at large xf in pA collisions 
are used to infer the existence of gluon depletion as the pro- 
jectile proton traverses the nucleus. The modification of the 
gluon distribution is studied by use of a convolution equation 
whose non-perturbative splitting function is determined phe- 
nomenologically. The depletion factor at xi = 0.8 is found to 
be about 25% at A= 100. 

PACS number: 25.75.Dw, 24.85. +p 



It is conventional in the study of J/^ production in 
heavy-ion collisions that the gluon distribution before the 
hard subprocess of cc production is assumed to be the 
same as in a free nucleon . The unconventional view 
that the gluon distribution can be modified in the nuclear 
medium due to depletion was suggested in In this 
paper we focus on p-A collisions and show that the data 
1^ on a{xp) can be used to infer that gluon depletion in 
the projectile proton is not negligible. 

Charmonium absorption in pA collisions has been stud- 
ied in §,|l without finding any satisfactory explanation 
for the xp dependence ofa{xF)- In Q the effect of energy 
loss of partons is considered, but that is only one aspect 
of gluon depletion. Here we pay particular attention to 
the evolution of the gluon distribution of the projectile 
as it traverses the nucleus. The approximate absence of 
dilepton suppression and the consequent implication that 
the quark distribution is nearly unaltered by the nuclear 
medium lead some to expect that the gluon distribution 
would be unaltered also. However, such a view is based 
on the validity of DGLAP evolution of the parton dis- 
tribution functions [1). We adopt the reasonable alter- 
native view that the evolution in a nucleus is different 
from that of pQCD at high Q^; indeed, we shall let the 
data guide us in determining the proper dynamics of the 
low-Q^ non-perturbative process. 

The Fermilab E866 experiment measured the J/^ sup- 
pression in p-A collisions at 800 GeV/c with a wide cov- 
erage of a;i? 1^. The result is given in terms of a{xp), 
which is defined by the formula 



a{xF) = 0.952(1 + O.G232:f - 0.397x|) 



(2) 



i?(xF, A) = aA{xF)/A(7NixF) = 



(1) 



where ctat^a is the cross section for J/ip production by a 
proton on a nucleon (N) or on a nucleus (A). In ||^ a 
paramctrization of ai^xp) for J/ip production is given: 



for —0.1 < xp < 0.9. It is our aim here to explore 
the implication of Eq. (|^) on the evolution of the gluon 
distribution. 

Since the semihard subprocess ofg-l-g— >c-|-cis com- 
mon for p-N and p-A collisions, they cancel in the ra- 
tio R{xf, A) so the xp dependence can come from three 
sources: (a) the ratio of the gluon distribution in the pro- 
jectile passing through a nucleus to that in a free proton, 
G{xf,A), (b) nuclear shadowing of gluons in the target, 
N{xf,A), and (c) hadronic absorption of the cc states 
after the semihard subprocess, H{xf,A). Putting them 
together, we have 



R{xF,A) ^ G{xF,A)N{xF,A)H{xF,A) 



(3) 



G(xf,^) and N{xf,A) are ignored in gj7[]. Since xp < 
0.25 in |q|, there is not much dependence on xp to be 
ascribed to H{xp,A), but in Q], where the full range 
of Xp is considered, H(xp,A) is forced to carry the en- 
tire a;i?-dependence by a fitting procedure, resulting in 
an unreasonably short octet lifetime. Our approach by 
including G{xf,A) and N{xfjA) in Eq.(||) is therefore 
complementary to the work of . 

The nuclear shadowing problem has been studied in de- 
tail by Eskola et al. [ [lO|jll ] , using the deep inelastic scat- 
tering data of nuclear targets at high Q^. On the basis of 
DGLAP evolution they can determine the parton distri- 
butions at any > 2.25 GeV^. The results are given in 
terms of numerical parametrizations (called EKS98 
of the ratio Nf-{x,Q'^) — fi/Ai^, Q'^)/ fiix, Q'^), where 
fi is the parton distribution of flavor i in the free pro- 
ton and fi/A is that in a proton of a nucleus A. We 
shall be interested in the ratio for the gluon distributions 
only at = 10 GeV^, corresponding to cc production, 
and denote it by N{x, A). From the numerical output of 
EKS98 we find that a simple formula can provide a good 
fit to within 2% error in the range 40 < ^ < 240 and 
0.01 <x < 0.12; it is 



N{x,A) = A>^'^''\ 



where 



f3{^{x)) = ^ (0.0284 + 0.0008 ^ - 0.0041 



(4) 



(5) 



with ^ = 3.912 -|- Inx. Thus the A dependence is minimal 
at ^ = 0, corresponding to a; = 0.02. 
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The variable x in Eq.(^ is the gluon momentum frac- 
tion in a nucleon in the nucleus, usually referred to as X2- 
Both xf in Eq.(P and X2 in Eq.(^ are to be converted 
to the xi variable for the projectile nucleon, using 



Xf = Xi - X2, 



xxX2=T = M']i^Is, (6) 



so that a part of Eq. (^) can be rewritten as 

i?(xj.,^)/iV(a;2,yl) = (7) 

In our approach we treat H(xf,A) as having negligible 
dependence on xf for all xp- Attempts [^,0 to find that 
dependence have failed and led to the suggestion of the 
existence of an unaccounted mechanism responsible for 
the enhanced suppression in R{xf, A) at large xp- In our 
view that mechanism is gluon depletion. Of course, if the 
Xp dependence of H{xf,A) were independently known, 
its incorporation in our analysis is straigtforward. For us 
here, we identify the xi dependence of G{xi,A) in Eq.(||) 
with that in Eq.(Q), which is completely known, and pro- 
ceed to the study of the phenomenological implication on 
gluon depletion. 

In the spirit of DGLAP evolution, even though the 
effect of a nuclear target on the projectile gluon distribu- 
tion is highly non-perturbative, we now propose an evo- 
lution equation on the gluon distribution g(x^ z), where z 
is the path length in a nucleus. For the change oi g{x,z), 
as the gluon traverses a distance dz in the nucleus, we 
write 



dz 



g{x,z) 



^ dx' X 
— 9ix',z)Q{ — ), 



X 



X 



(8) 



where Q{x/x') describes the gain and loss of gluons in 
dz, but unlike the splitting function in pQCD, it cannot 
be calculated in perturbation theory. Equation (||) is 
similar to the nucleonic evolution equation proposed in 
[ p2[ , except that this is now at the parton level. Instead 
of guessing the form of Q{x/x'), which is unknown, we 
shall use Eq.(|^) to determine it phenomenologically. 
To that end, we first define the moments of g{x, z) by 



(9) 



gn{z) = / dxa;" g{x,z). 
Jo 

Taking the moments of Eq. (0) then yields 



dgn{z)/dz = g„(z) (3„, (10) 
where Qn — dyy"^^"^ Qiu)- It then follows that 

<?n(^)=5n(0)e^«", (11) 

whose exponential form suggests Qn < for the phys- 
ical process of depletion. The gluon depletion function 
D{y, z) is defined by 



g{x,z) = I ^g{x',0)D{^,z), 
x' x' 



(12) 



where g{x' , 0) is the gluon distribution in a free nucleon. 
From Eq. (|l|) we have g„(z) = g„(0) D^iz), where Dn{z) 
is the moment oiD{y, z) . Comparison with Eq. (O) gives 



Dn{z) 



(13) 



To relate this result to R{xp,A), we first note 
that G{xp,A) in Eq.(||) is, by definition, G{xp,A) — 
g{xi, A) / g{xi,Q), where xp is expressed in terms of xi. 
It then follows from Eq.(|) that 

J{xi,A) = g(a;i, 0) R{xp{xi), A)lN{x2{xi),A) 

^g(xi,A)H{A). (14) 

In relating A to the average path length L of the pro- 
jectile p through the nucleus, we use L = 3i?A/2 — 
l.SA^/^bn. We then set z = L/2 for the average dis- 
tance traversed at the point of cc production. Thus 
when referring to the last expression of Eq.(|lJ), we write 
J{xi,A) = g{xi,z{A)) H{z{A)), where g{xi,z) is to be 
identified with that in Eq.(|lJ). Note that the A depen- 
dence of the middle expression in Eq.(^) is, on account 
of Eq.(^, in terms of \nA, whereas that of the last ex- 
pression is in terms of z, or A^^^. Since it is known that 
\nA « ^1/3 for 60<A< 240, we shall consider the con- 



sequences of Eq.(|14|) only for A in that range. We suggest 
that a revised form of presenting the data, different that 
in Eq.(0), should be tried in the future. 

Taking the moments of J{xi,A), we get using Eq.(ll) 



lnJ„(yl) - lng„(0) = zQ„ + lnH{z). 



(15) 



To determine Qm it is necessary to use as an input 
the gluon distribution g{xi,0) in a free proton at = 
lOGeV^. We adopt the simple canonical form 



gixi,0) ^ goil - xif 



(16) 



where the constant go is cancelled in Eq.(|15D due to the 
definition of J{xi,A). In our calculation we set go = 1. 
Indeed, the accuracy of g{xi,0) is unimportant, since it 
enters Eqs.(|l|) and ^ in ways that render the result in- 
sensitive to its precise form. On the basis of Eqs.(0) and 
(|l|), J{xi,A) is therefore known. The LHS of Eq.(|l|) 
can then be computed except for a caveat. To calcu- 
late the moments of J{xi,A), it is necessary to compute 
dxi ^(a^i. A). However, xi cannot be less than t 
in order to keep X2 < I [see Eq.(||)]. Furthermore, Eq.(j7) 
does not provide reliable information on J{xi , A) at small 
xi, since the parametrizations of a{xp) and P{x2) are for 
the variables in ranges that exclude the xi — s- r limit. 
Fortunately, that part of the integration in xi can be 
suppressed by considering n > 3. The part of the inte- 
gration in the interval < xi < r amounts to only about 
2% contribution even at n = 2 (if naive extrapolation is 
used), so its inaccuracy will be neglected. Physically, it 
is the data at high x p that we emphasize in our analysis, 
and that corresponds to the high-n moments of J{xi, A). 
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For convenience, let us denote the LHS of Eq.(]l5|) by 
Kn{z), i.e., Kn{z) EE ln[J„(z(A))/5„(0)]. For sample 
cases of j4 = 100 and 200, they are shown as discrete 
points in Fig. 1 for 3 < n < 20. Instead of perform- 
ing an inverse Mellin transform on Kn{z), our procedure 
is to fit Kn{z) by a simple formula that can yield Q{y) 
by inspection. The fitted curves shown by the solid and 
dashed lines in Fig. 1 are obtained by use of the formula 



-fco H 

n n 



k-2 



1 



(17) 



Using ki and k[ to denote the values for the cases A — 100 
and 200, respectively, we have 

A:o = 1.592, ki = 23.42, fca = 97.66, fcg = 89.17 
/co = 1.831, k[ = 27.43, k'^ ^ 113.97, fcg = 103.80. 



Because of Eq.(|l5|), the n dependence of Kn prescribes 

(18) 



the n dependence of Qn- Let us therefore write 
91 92 , 93 



-90 
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Since Eq.(p5|) is to be used only for A > 60, we evaluate 
it at ^ = 100 and 200, and take the difference. Denoting 
z by zi and Z2, respectively, for the two A values, and 
with Aki ^ k[ — ki, Az = Z2 — zi, we have 



Afco ~ 9o — In 



H{Z2) 

H{z,y 



Ah^q.Az, (i^O). (19) 



For the hadron absorption factor H{z) we write it in 
the canonical exponential form p3[ |, H{z) — exp(— pcrz), 
where p-^ = (4/3)7r(1.2)3 fm^, z = 0.9A^/^fm, and a is 



the absorption cross section. Putting these in Eq.(19), 
we get (with Az — 1.086 fm) 

qo + pa = 0.22, qi = 3.68, q2 = 15.01, 93 = 13.47 (20) 

in units of fm~^. 

There is a reason why go SLiid pa enter Eq.(pO|) as a 
sum. To appreciate the physics involved, we first note 
that Eq.(|l^) imphes directly 



Qiv) -90 <5(1 -y) + qiy-q2 z/ + 93 V^- 



(21) 



The first and third terms on the RHS above are the loss 
terms (i.e., gluon depletion), while the second and last 
terms represent gain (i.e., gluon regeneration). If Q{y) 
consisted of only the first term, then using it in Eq.(|^) 
would give dg(x, z)/dz — —qo g{x, z), whose solution is of 
the same exponential form as that of absorption. With 
both depletion and absorption present, the exponents 
lead to a sum, as in Eq.(po[). Our Q{y) is, however, more 
complicated. The —52 y^ term gives rise to depletion that 
depends on the shape of g{x,z), while the qiy + qsy^ 
terms generate new gluons at x from all the gluons at 
x' > X. 

Since Qn decreases monotonically with n, we require 
Qs < 0, and exclude Q2 from this consideration because 



of its inaccuracy discussed earlier. Combining Eqs.(18) 
and (20), we get pa < qo + po' — 9i/3 + 92/4 - 93/5 = 0.05 
fm"^ We thus set qo = 0.17 fm"^ 

Since it is not easy to see directly from Qn or Q{y) the 
magnitude of the effect of gluon depletion and regenera- 
tion, we can calculate g{xi, z), not from Eq.(|l^), but by 
fitting the calculated gn{z) in Eq.(ll), using the formula 



9n{z) = ^ ai(z) B{n -1,5-1- i), 



(22) 



where B{a, b) is the beta function. Then the result yields 
directly 



g{xi,z) = ^ ai{z) (1 - xi) 



4+i 



(23) 



For A — 100 (200), i.e., z = zi (Z2), we have ai = 
0.58 (0.485), as = 0.92(1. 118), and 03 = -0.47 (-0.56) 
for go = 1 in Eq.(^. The result for G{xi,z) = 
g{xi, z)/ g{xi,0) is then 

G{xi,z) = ai(z) -I- a2(z)(l - xi) + a3(z)(l - xif, (24) 

which is shown in Fig. 2 for two values of A. It is now 
evident that gluon depletion suppresses the gluon dis- 
tribution at medium and high xi, but the unavoidable 
gluon regeneration enhances the distribution at low xi. 
The cross-over occurs at xi « 0.28. 

Let us now exhibit our result for a{xF), which is shown 
in Fig. 3. The line is obtained by use of Eq.(|4|) in Eq.(||) 
and a = 6.5 mb in H{z). Only one line is shown for 
both A = 100 and 200, their difference being negligible 
in the plot. Since our method of using the moments can- 
not be extended to n = 2 due to the problems mentioned 
after Eq.([T^), there is some inaccuracy inherent in our 
analysis. Thus the fit cannot be expected to be perfect. 
Our model can reproduce the general trend of the xp 
dependence, but not the detail structure, for which more 
terms in Eqs.(p7|) and ( |l8| ) would be needed. The over- 
all suppression is achieved by use of a phenomenological 
value of a, rather than the bound based on the technical 
assumption of Q3 < 0. 

Our analysis has been based on the assumption that 
H(z,A) is independent of a;^- If and when that xp de- 
pendence can be determined independently, the effect can 
easily be incorporated in our analysis to modify our nu- 
merical result . Since that dependence is not likely to be 
strong |^-|^, the modification would be minor. Our study 
shows that the J ftp suppression observed at large xp 
pA collisions strongly suggests the presence of gluon 
depletion in the beam proton at high xi . The significance 
of this finding goes beyond the J/ suppression problem 
itself, since it would revise the conventional thinking con- 
cerning the role of partons in nuclear collisions. 

Since the gluon distribution is enhanced for xi < 0.28, 
the J/ip suppression observed in the xp ^ region in the 
heavy-ion collisions at CERN-SPS cannot be due to the 
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gluon depletion effect. The same would be true at RHIC. 
However, we expect a significant increase in suppression 
at large xp due to gluon depletion, not to color decon- 
fincmcnt. We further speculate at this point that the 
gluon enhancement at low x may be responsible, at least 
in part, for the strangeness and dilepton enhancement 
already observed in heavy-ion collisions. 

This work was supported, in part, by the U.S. -Slovakia 
Science and Technology Program, National Science Foun- 
dation under Grant No. INT-9319091 and by the U. 
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FIG. 1. K„. The curves are fitted results using Eq.(17). 



FIG. 2. G{xi,z) showing the effects of gluon depletion. 

FIG. 3. a{xF) vs xf- The solid line is our result compared 
to the data from [3]. 
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